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Abstract 

In this paper, we consider the energy decay of a damped hyperbolic system of wave- 
wave type which is coupled through the velocities. We are interested in the asymptotic 
properties of the solutions of this system in the case of indirect nonlinear damping, i.e. 
when only one equation is directly damped by a nonlinear damping. We prove that the 
total energy of the whole system decays as fast as the damped single equation. Moreover, 
we give a one-step general explicit decay formula for arbitrary nonlinearity. Our results 
shows that the damping properties are fully transferred from the damped equation to 
the undamped one by the coupling in velocities, different from the case of couplings 
through displacements as shown in [DlllIIl for the linear damping case, and in [7] for 
the nonlinear damping case. The proofs of our results are based on multiplier techniques, 
weighted nonlinear integral inequalities and the optimal-weight convexity method of [Si- 

Keywords Energy decay, nonlinear damping, hyperbolic systems, wave equation, plate equa¬ 
tion, weighted nonlinear integral inequality, optimal-weight convexity method. 

1 Introduction 

Let n be a bounded subset of R"' with a smooth boundary denoted by E. We consider 
the following wave system 

/ 

tt" — Au-b a(x)u'-b/o(x, u') = 0 in nx(0,-boo), 
v" — Av — a{x)u' = 0 in n X (0, -boo), 

u = u = 0 on r X (0, -boo), 

(n, u')(0) = (u°, u^), (u,u')(0) = (u°,u^) in O . 
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To a strong solution of this system, we associate the energy defined by 

E{t) = l [ (|nf + |Vw|2 + \v'\^ + \Vv\‘^)dx . (1.2) 

Jn 

One can show that the energy of the strong solutions of this system satisfies 

E'{t) = — / v!p{x,u')dx , (1.3) 

Jn 

In the sequel, the assumption on p will ensure that /?(•, s)s ^ 0 a.e. in and for all s £ M, 
so that the energies of the strong solutions satisfy E'{t) < 0. Hence the nonlinear term in 
our coupled system is indeed a damping term, so that one expects the energies to decay to 0 
at infinity. 

Let us recall that for the scalar damped wave equation, that is for 

u” — Au + p{x, u') = 0 in H X (0, +oo), n = 0 on T x (0, +oo). 

It is well known that when the damping term is linear, i.e. when p{.,u') = a[-)u', where 
a ^ 0 a.e. on H, the energy of the solution decays exponentially under some some geometric 
conditions on the support of a (see [13 EH ). When the damping term p{-,-) is nonlinear 
with respect to the second variable, [5] and m, m give, under some suitable geometric 
conditions, a one-step explicit energy decay formula in terms of the behavior of the nonlinear 
feedback close to the origin. These results rely on a general weighted nonlinear integral 
inequality together with an optimal-weight convexity method developed in [5] ■ If no geometric 
assumptions on the damping region are made, the decay is known to be of logarithmic type 
for a linear damping (see e.g. mm)- 

Let us go back to the above wave-wave coupled system m- At this stage, four main 
features characterize this system: 

• only the first equation is damped 

• the damping p may have an arbitrary growth around 0 with respect to the second 
variable 

• the coupling coefficient a may vanish in some parts of H 

• the coupling is acting through the velocities. 

Let us now comment on these features. 

• The fact that only one equation of the coupled system is damped refers to the so- 
called class of ” indirect” stabilization problems initiated and studied in mmm and further 
studied in ig la 113]. Indeed, when dealing with coupled systems, it may be impossible 
or too expensive to damp each equation. Such an example is provided for instance by the 
Timoshenko system [7iEgEn|. More generally, coupled system involving some undamped 
equations, are said to be indirectly damped. From the point of view of applications in control 
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theory, a challenging question is to determine whether the single feedback is sufficient to 
guarantee that the energy of the full system decays to 0 at infinity and to determine at which 
rate. In this latter case, the lack of feedback on the second equation is compensated by the 
coupling effects. 

• The case of general damping feedbacks, that is with arbitrary growth close to 0, has 
received a lot of attention since more than a decade. The first result in this direction has 
been derived in [22], however no general simple explicit formula was provided except for 
linearly or polynomially growing dampings close to 0. Such first examples of explicit general 
formula are given in m (see also [26] ), but this formula does not allow to recover in a single 
step the expected quasi-optimal energy decay rate in the polynomial case (or for polynomial- 
logarithmic growth). As far as we know, the first result giving a general one-step quasi-optimal 
semi-explicit formula is given in [5]. A further analysis based on a suitable and original 
classification of the feedback growth has been introduced in [8]. This classification gives 
a very simple one-step explicit energy decay formula for general feedbacks growth, provided 
that this growth is not close to a linear behavior. A more complex semi-explicit formula holds 
in the general case including feedbacks with a growth close to a linear one around 0. By a 
one-step formula, we mean here, a formula which gives a decay rate depending explicitly on 
the feedback in a simple explicit way. In particular, this formula does not require further steps 
as in most of the existing literature to lead to explicit expressions. Moreover the optimality of 
the formula is proved for the corresponding finite dimensional systems or for semi-discretized 
scalar wave or plate equations, whereas optimality results are proved for some examples in 
the infinite dimensional case in [5], using results of [311 [32]. For previous results and also 
results for semilinear damped equations, one can see [niEiiEaisiiM]. 

• When the coupling coefficient is bounded below by a strictly positive constant, the 
coupling is active in the whole domain fl, so that the equations are coupled in the whole 
domain. If a is nonnegative on 12 but is allowed to vanish on a subset of 12, the equations 
are ’’uncoupled” in I2\supp{a}, so that in this region the second equation is decoupled from 
the first equation and is undamped. Such cases are harder to handle. A hrst study in this 
framework, but for linearly damped wave-type equations coupled in displacements is given 
in [T3] . 

• When indirect damping occurs through displacements, that is for systems coupled in 
displacements, as for instance 

u" — Au -|- a{x)v a{x)u' = 0 in (0, oo) x 12, 
v" — Av -\- a{x)u = 0 in (0, oo) x 12, 
u = u = 0 on (0, oo) X r, 

(u(0), u'(0)) = (u°, u^), (u(0), ?;'(0)) = (u°, u^), in 12, 

it has been shown in [2] that even for constant coefficients a, the energy of this linearly 
damped system never decays exponentially, but decays only polynomially with a decay rate 
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depending on the smoothness of the initial data and a general lemma announced in [T] (see 
mm for a proof). These results are based on the method of higher order energies initiated 
in [1] and developped in [2] and [3] for the indirect boundary damping cases in an abstract 
setting and applied to various examples. The result of [2] has been generalized to coefficients 
a that vanishes on a subset of in [T3], under certain assumptions on the supports of a 
and a (roughly speaking they are both supposed to satisfy the Geometric Control Condition 
(GCC) of [IS]). Further results on coupled models with distributed dampings but satisfying 
hybrid boundary conditions have been obtained in [T2|. Sharper results have been obtained 
through interpolation techniques extending the first results of [T]. Let us further mention 
that in m, the author shows that the energy decays logarithmically under the assumption 
supp{a} n supp{a} / 0. Hence when the coupling acts through displacements, indirect 
stabilization occurs but in a weaker form than the one of the corresponding scalar case, since 
exponential stabilization does not hold even for a linearly damped case. 

The goal of this paper is to generalize the quasi-optimal energy decay formula given for 
scalar wave-type systems in |H] (see also 0), to the case of coupled systems in velocities, 
under the above four features. More precisely, we prove, under some geometric conditions 
on the localized damping domain and the localized coupling domain, that the energy of 
these kinds of system decays as fast as that of the corresponding scalar nonlinearly damped 
equation. Hence, the coupling through velocities allows a full transmission of the damping 
effects, quite different from the coupling through the displacements. 

The optimality of the above estimates has been proved for finite dimensional equations, 
including the semi-discretized wave equations in |8] . In the inhnite dimensional setting, lower 
energy estimates or optimality are open questions. Optimality has been only proved in the 
particular case of one-dimensional wave equation with boundary damping (see ISlEIllSI). 
Lower energy estimates have been established in |8l[9lll0| for scalar one-dimensional wave 
equations, scalar Petrowsky equations in two-dimensions and Timoshenko systems. We use 
the comparison method developed in |10j to extend these results to one-dimensional wave 
systems coupled by velocities. 

Remark 1.1. The method presented here is general and can easily be adapted to handle 
corresponding coupled systems of plate equations, elasticity models, and more complex ex¬ 
amples in the spirit of the general approach given in [5]. Our aim through this paper is to 
give a general methodology on a concrete PDF example to show that if the damping effects 
are suitably transferred through the coupling operators, then indirect stabilization can pro¬ 
duce damping mechanisms of the quality of a direct damping for the corresponding scalar 
equation. 

Note that one can also present, with no additional mathematical originality and no gain 
with respect to applications, all our results by means of a more ’’Lyapunov” presentation. In 
this case, it is sufficient not to integrate over the time interval and to handle terms of the 
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form ^ for instance, and multiply afterwards by a weight function, which can be 

a weaker (and less good weight function) than in the original method introduced for the first 
time in [5] and announced in [3]. 


The rest of this paper is organized as follows: In Section [2l we give some basic preliminar¬ 
ies, assumptions and notations. The main results, including Theorems 13 .1 1 on well-posedness 
of (|l.ll) , Theorem 13.21 on energy decay for polynomially growing damping case. Theorem 
[331 on energy decay for general nonlinear damping case. Theorem 13.41 on lower energy es¬ 
timates for one dimensional system (1321), are presented in Section [3l Explicit decay rates 
corresponding to some typical dampings are also provided in Subsection 13.41 As the main 
tool of deriving the quasi-optimal one-step explicit energy decay formula, the optimal-weight 
convexity method together with general weighted nonlinear integral inequalities are intro¬ 
duced in Section 01 The proof of the main results, as well as the decay rates of Example 
13.1113.41 are given in Section [5l Finally, Section [6l is used to prove Lemma 15.11 on weighted 
energy estimates for a single non homogeneous wave equation by the multiplier method. 


2 Preliminaries, assumptions and notations 

Notations 

For the simplicity of statement, we denote in the whole paper L^(0) by L^, Hq{Q) 
by Hq, by Moreover, we say that the initial data are in the energy space 

whenever E Hq x and G Hq x and the initial data are smooth if 

G n H^) X and G n H^) x 

Geometric conditions 

As already mentioned in the introduction, stabilization results for wave-like systems 
require geometric conditions on the region where the feedback is active. In the sequel, 
we shall consider the so-called Piecewise Multipliers Geometric Condition (denoted by 
PMGC, in short): 

Definition 2.1 (PMGC). We say that a subset cu C 11 satisfies the PMGC, if there 
exist subsets Qj C H having Lipschitz boundaries and points xj G M'^, j = 1, • • • , J, 
such that n = 0 for i ^ j and uj contains a neighborhood in H of the set 
Uj=i 7 j(xj) U (11 \ where ^j{xj) = {x G clHj : {x — Xj) ■ Vj{x) ^ 0} and Vj is 

the outward unit normal to (911j. 

Assumptions on the coupling coefficient 
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We assume that the coupling function a € (7(0) satisfies 

3 a+ > 0, a_ > 0 such that 

< Q!+ ^ a{x) ^0, Vx G O, (HC) 

a{x) ^ a- > 0, Vx G Wc C O, 
where and ujc is an open subset of O with positive measure. 


Assumptions on the feedback 


We consider feedbacks p with an arbitrary growth close to 0. However, to give the reader 
a better insight of the scope and challenge of one-step explicit general quasi-optimal 
energy decay formulas, we hrst provide the result and proof for polynomially growing 
feedbacks, for which the proofs are easier and then the general result for arbitrary 
growing feedbacks. Hence, we detail below the two sets of assumptions: the one for the 
polynomial case, then those for the general case. 

The assumptions in the case of polynomially growing feedbacks is as follows 

p e C(fl X R), p{x, 0) = 0 Vx G H , 

s !-)■ p{x, s) is nondecreasing V x G H, 

3 c > 0 and p ^ 1 ,3 a G C{Q) such that 

a(x)|s| ^ \p{x, s)| ^ co(x)|s|, Vx G H , |s| ^ 1, 

1 

a(x)|s|P ^ \p{x, s)| ^ ca(x)|s| p, VxGH,|s|^ 1, where 

a^0onH,3a_ >0 such that a(x) ^ a_, Vx G C H, 
where is an open subset of H with positive measure. 


The assumptions in the case of arbitrary growing feedbacks is as follows 


/9 G (7(H X M), p{x, 0) = 0 Vx G H , 
s !->• p{x, s) is nondecreasing V x G H, 

3 c > 0,3 a G C{Q) and 3 g( G such that 

< a(x)|s| ^ |/o(x, s)l ^ ca(x)|s|, Vx G H , |s| ^ 1, 

a(x)( 7 (|s|) ^ |/9(x, s)| ^ ca(x)( 7 “^(|s|), VxGH,|s|^ 1, where 

a^0onH,3a_>0 such that a(x) ^ a_, Vx G C H, 

(7 is a strictly increasing and odd function. 


Remark 2.1. Thanks to the hypotheses (HFp) or (HF„), we have 


p{x, s)s ^0, VxGH,VsGM, 


(HF,) 


( 2 . 1 ) 


which ensures that the energy of the solutions of the above wave system is nonincreasing. 

Remark 2.2. Note that we can infer from [HFg) that for very e G (0,1), there exists 
constants ci > 0, C 2 > 0 such that 

I cia(x)|s| ^ |p(x,s)| ^ C2a(x)|s|, VxGH,|s|^e, 

[ cia(x)5r(|s|) ^ \p{x, s)| ^ C 2 a{x)g~^{\s\), Vx G H , |s| ^ e . 
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Convexity assumptions on the feedback and some definitions 


Assume that (HF^) holds. Then, following mm we assume that the function H defined 
by 

H{x) = y/xg{y/x), (2.3) 


is strictly convex in a right neighborhood of 0, i.e. on [0, Tq] for some sufficiently small 
ro € (0,1]. We define the function i? on M by H{x) = H{x) for every x € [0,rQ] and 
by H{x) = +00 otherwise, and we define the function L on [0, +oo) by 


L[y) 


' H*{y) 

' y 

0 


if y > 0, 


if y = 0, 


(2.4) 


where H* is the convex conjugate function of Ff, defined by 


#*(y) = sup{xy - H{x)}. 

xGM. 

By construction, the function L : [0,+oo) ^ [0, Tq) is continuous, one-to-one, onto and 
increasing, moreover it is easy to check that 


0 < < ri 


(2.5) 


holds (see mm for a complete proof). We also define the function Ah on (OjTq] by 


Ah{x) 


H{x) 
xH'{x )' 


( 2 . 6 ) 


Remark 2.3. The function Ah has been introduced for the first time by the first author 
in [8]. It is an essential tool to classify the feedback growths around 0 and to simplify 
the decay estimate formula given in -without loosing optimality properties- for the 
feedbacks having a growth around 0 which is not close to a linear one (as explained 
below). 


Remark 2.4. Note that due to our convexity assumptions, we have 


Ah{[0,1]) C [0,1]. 


When g{x) = x (linear feedback case), Ah{-) = 1. If we now set for instance, g{x) = 
x(ln(l/x))“'? in (0,e] with q > 0 and e > 0, then hmsupA/f(x) = 1. Many other 

a;->0+ 

examples of feedbacks such that limsupA//(x) = 1 can be given, they characterize a 

3;^0+ 

growth which is close to a linear one. This leads to the following definition: 


Definition 2.2. We say that a feedback p satisfying (HFg) has a growth close to a 
linear one in a neighborhood of 0, if it is such that the function H defined by (12.3|) 
satisfies limsupA//(x) = 1. Otherwise, one says that the feedback p is away from a 

a:^0+ 
linear growth. 
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On the opposite side, for functions g which converge very fast to 0 as x goes to 0, 
such as for instance g{x) = for x G (0,e] (and many other examples), one has 

limsup A//(x) = 0. 

3;^0+ 

For polynomially growing feedbacks, e.g. when g{x) = with p > 1, we have 
Ah{-) = For feedbacks such as g{x) = x^’(lnl/x)'^ with p > l,g > 0, we still 

have limsup A/i-(x) = 

x^0+ 

We will see later on, in Theorem ld.dl that the case of feedbacks close to a linear behavior 
as X goes to 0 has to be distinguished from the other cases. 


Finally, we define, for x > l/H'{rQ), 


(2.7) 


Remark 2.5. Note that when g'{0) 0, g has a linear growth close to 0. Therefore, 

this case is similar to the linear case which is already well-known. We thus focus in the 
sequel on the cases where 5 ^( 0 ) = 0. 


3 Main results 


3.1 Well-posedness 

We set n = X and set U = {u,p,v,q). We equip T-L with the scalar product 


{U ,U) = I [Vu-Vu + pp + Vv-Vv + qqjdx 


Theorem 3.1. Assume (HF^) and that a G Then for all initial data in energy 

space, there exists a unique solution {u,v) G C([0,+oo); (Rq)^) nC^([0,+oo); (L^)^) of (|l.ll) . 
Moreover for any smooth initial data, the solution satisfies {u, v) G L°°([0, +oo); (R^nilQ)^)n 
VF^’°°([0,+oo); (Rq)^) n kF^’°°([0,+oo); (L^)^). Moreover, in this latter case the energy of 
order one, defined by 


IS non increasing, i.e 


= \ [ {\utt\‘^+ \Vut\'^+ \vtt\‘^+ \Vvtf)dx , 

^ Jn 

■} 

Ei{t) < F'i(O). 


(3.1) 

(3.2) 


3.2 One-step quasi-optimal energy decay rate for the wave-wave system 
First case: polynomiaily growing dampings close 0 

For the sake of clarity, we first provide the results in the case of a polynomially growing 
feedback s ^ p{-,s). 
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Theorem 3.2. Assume that p > 1, (HFp) and dHCj hold. Assume also that Ud and lOc 
satisfy the PMGC. Then there exists a* > 0 such that for any G (0,a*] and any non 
vanishing initial data in the energy space, the total energy of dLU) defined in decays as 

E{t) ^ t^, Vi € [Te(o) , +oo) , (3.3) 

where Ce{o)-,Te[o) > 0 are constants depending on E{{)). 

Remark 3.1. If p = 1, we can follow the proof of Theorem 13.21 and obtain, by using Lemma 
14.21 instead of Lemma l4.ll the exponential stability for system (11.11) 

E{t) ^ CE{id) e-^\ Vi G [0, +oo), (3.4) 

where (7 > 0, k > 0 are constants independent of the initial data. 

Second case: arbitrarily growing dampings close to the origin 

If we consider a more general nonlinear damping p, we provide below a quasi-optimal one- 
step explicit energy decay formula following the optimal-weight convexity method together 
with general weighted nonlinear integral inequalities developed in UM- 

Theorem 3.3. Assume (HFg) and (jHCp hold. Assume that the function El defined by (12.3p is 
strictly convex in [OjTq] for some sufficiently small ro G (0,1] and satisfies H{0) = H'{0) = 0. 
We define the maps L and Ah respectively as in (Eap and (1^ . Assume also that ujd and 
ujc satisfy the PMGC. Then there exists a* > 0 such that for any a_|_ G (0,a*] and any non 
vanishing initial data in the energy space, the total energy of dLU) defined in dLl decays as 


E(t) < , V t 

^0 V M 1 


> 


M 


H'{r, 


2 ^ ’ 


(3.5) 


where I3 e{o) defined by (|5.30l) . M is defined by (I5.34p and independent of E{0). Further¬ 


more, i/limsup A/f(a;) < 1, then E satisfies the following simplified decay rate 

x^0+ 

1 / kM \ 


E{t)<2/3E(o)(H' 


t 




(3.6) 


for t sufficiently large, and where k> D is a constant independent of EifS). 

Remark 3.2. Note that when g{x) = wih p > 1, then H{x) = so that 

limsupAj:/(x) = —^ < 1. Then the formula ()3.6I) gives back the energy decay rate of 

i-2/(p-i) given in Theorem 13.21 


Remark 3.3. The smallness of a can be reduced if we assume additionally supp{a} C iWd 
in Theorem 13.31 Actually, we can choose <5 = 1 in the proof of (IS.lUp . and in that case, the 
last term of (Ib.lOp can be replaced by 

Cl f 4>{t) f \u'\‘^dxdt 

J S J 

Then (|5.3p follows easily by choosing ei > 0 sufficiently small in (15.101) . The proof afterwards 
is the same. 
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3.3 Lower energy estimates 


The optimality of the above estimates are open questions. Here we use the comparison 
method developed in [TO] to establish the lower estimate of the energy of the one-dimensional 
coupled wave system. Actually, we consider the following wave-wave system in O = (0,1) C M 

and p{x, s) = a{x)g{s) for all x G 0 and all s G M. 

/ 

utt-Uxx + ct{x)vt + a{x)g{ut) = 0, in Qx{0,+oo), 

vtt - Vxx - a{x)ut = 0 , in n X (0, -hoc), 

u{t, 0) = u{t, 1) = v{t, 0) = v{t, 1) = 0 , for t G (0, -|-oo), 

{u,ut){0,x) = {u^{x),u^{x)), {v,vt){0,x) = {v^{x),v^{x)), for xGtl. 

•. 

We define H by (12.3p . Ah by (|2.6I) and consider the following assumptions 


3ro > 0 such that the function H : [0, rg] >->• M defined by (12. 3p 
is strictly convex on [0, rg] and 

either 0 < liminf Aff(x) < limsup Aj;f(a;) < 1 ('HF;i 

or there exist g > 0 and zi G (0, zq] with (|5.39l) such that 
° 'm Hmsup Aj^(x) < 1. 

Theorem 3.4. Assume that (HF^), ( |HF/ ), and ()HCp hold. Assume also that ojd and lOc 
satisfy the PMGC, and that a+ G (0, a*] where a* > 0 is as in Theorem \3.3[ Then for all 
non vanishing smooth initial data, there exist Tq > 0 and Ti > 0 such that the energy of (13771) 
satisfies the lower estimate 


E{t)> 


1 




H') 


-1 


1 


-t — Tq 

where js = 4y^Ei(0), Cg is as in Lemma El 


V t > Ti + To 


(3.8) 


3.4 Some examples of decay rates 

For the sake of completeness, we give some significative examples (taken from jS] [8]) of 
feedback growths together with the resulting energy decay rate when applying our results. In 
the sequel, Ce{o) > 0 stands for a constant depending on E{Q), while C'^j(g) > 0 is a constant 
depending on Ti(0). 

Example 3.1. [The Polynomial Case] Let g {x) = on (0, rg] with p > 1. Then the energy 
of (11.11) satisfies the estimate 

E{t) < CE(o)t~~, 

for t sufficiently large and for all non vanishing initial data in the energy space. Moreover, 
the energy of (|3.7p satisfies the estimate 

E{t) ^ 
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for t sufficiently large and for all non vanishing smooth initial data. 

_ 1 

Example 3.2. [Exponential growth of the feedback] Let g{x) = e ^ on (0,ro]. Then the 
energy of (HTD decays as 

E{t)<CEi^o) (ln(t))-', 

for t sufficiently large and for all non vanishing initial data in the energy space. Moreover, 
the energy of (|3.7p satisfies the estimate 

for t sufficiently large and for all non vanishing smooth initial data. 

Example 3.3. [Polynomial-logarithmic growth] Let g (x) = (In (l/x))*^ on (0, ro] with 
p > 2 and g > 1. Then the energy of a decays as 

E{t) < 

for t sufficiently large and for all non vanishing initial data in the energy space. Moreover 
the energy of (|3.7p satishes the estimate 

for t sufficiently large and for all non vanishing smooth initial data. 

Example 3.4. [Faster than Polynomials, Less than Exponential] Let g (x) = on 

(0, ro] with p > 2. Then the energy of (11.11) decays as 

for t sufficiently large and for all non vanishing initial data in the energy space. Moreover 
the energy of (|3.7p satisfies the estimate 

E{t) > 

for t sufficiently large and for all non vanishing smooth initial data. 

Remark 3.4. For these above four examples, one can show that lim Ah (x) < 1. Moreover, 

3:^0+ 

it is proved in [H], that the above decay rates are optimal in the finite dimensional case. 

4 Weighted nonlinear integral inequalities and decay rates 

Definition 4.1. We say that a nonnegative function E defined on [0, -|-oo) satisfies a weighted 
nonlinear integral inequality if there exists a nonnegative function w defined on [0, r/) with 
0 < T] < -|-oo and a constant M > 0 such that E([0, +oo)) C [0, g) and 

w{E{s))E{s)ds < ME{t) V t ^ 0 . (4.1) 
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Definition 4.2. We say that such a weight function w is of ’’polynomial” type, if there exists 
a > 0 such that 

w{s) = V s ^ 0, with 7] = + 00 . 

It is well-known that when E is a nonnegative, nonincreasing absolutely continuous function 
satisfying (HU) with a polynomial weight function, then E satisfies an optimal decay rate at 
infinity, proved in Ell Theorem 9.1] (see also herein for other references), that we recall in 
the next subsection. 


4.1 Polynomial weights 

Lemma 4.1. \21i Theorem 9.1] Assume that E : [0,-|-oo) [0,-too) is a non-increasing 

function and that there are two constants a > 0 and T > 0 such that 

OO 

E'^+\s)ds < TE{0)‘^E{t), yt G [0, +oo). 

Then we have ^ 

< ^(0) (1^) VtG[r,+oo). 

Remark 4.1. Note that when applied to decay estimates for dissipative systems, the above 
Lemma has to be used with a constant T which blows up as E[Q) = 0, so that the minimal 
time for which the above decay estimate is valid, blows up as £'(0) = 0. This result can be 
reformulated as below, to give an estimate which is valid for E{tf) > 0 as well as for E[0) = 0 
and for any t ^ 0 as explained below. 


Corollary 4.1. Assume that E : [0, -too) e- >■ [0, -too) is a non-increasing function and that 
there are two constants a > 0 and M > 0 such that 


E‘^^\s)ds < ME{t), VtG[0,+oo). 


Then we have 


E{t) < E{0) min V t ^ 0. 

\\M + aE(0)'^tJ J 


.M + aE{0)'^t. 

In particular for E{0) > 0, we deduce that 


E{t) < E{0) 


M(a-tl) \i/“ 
M + aE{{)Yt) 


V t ^ ME{0)-^/^ . 


Lemma 4.2. 1211 Theorem 8.1] Assume that E : [0,-too) [0,-too) is a non-increasing 

function and that there is a constant T > 0 such that 



E{s)ds < TE{t), 


Vt G [0, -too). 


Then we have 

E{t) < E{0)e^-^/^, yt G [T, +oo). 
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4.2 General weights 


For general weight functions, semi-explicit optimal decay rates have been derived for the 
first time in [5], and later on a simplihed form of the rates in [8]. 

Let r] > 0 and M > 0 be hxed and w he a strictly increasing function from [0, r]) onto 
[0,-|-oo). For any r G (0, ?]), we define a function Kr from (0,r] on [0,-|-oo) by; 

dy 


/ r 


yw{y) ’ 


(4.2) 


and a function ijjr which is a strictly increasing onto function defined from -|-oo) on 

'tpr{z) = Z + Kr{w~^{-)) > Z , V Z > —4 . (4.3) 

z wyr) 

We can now formulate our weighted integral inequality: 

Theorem 4.1. |3l Theorem 2.1] We assume that E is a nonincreasing, absolutely continuous 
function from [0,-|-oo) on [0,r]), satisfying the inequality 


f 


w{E{t))E{t)dt <ME{S), V0<5<r. 


Then E satisfies the following estimate: 

E{t) < w-^ 

where r > 0 is such that 

n+oo 


V’r {m) 


V t > 


M 

w{r) ’ 


(4.4) 


(4.5) 


^ + CX> 

— / E{t)w{E{t)) dr < r < T]. 

J 0 


M 

In particular, we have ^ lim E{t) = 0 with the decay rate given by ra- 

Theorem 4.2. Theorem 2.3] Let H be a strictly convex function on [0, rg] such the 
H{0) = H'{0) = 0. IFe define L and Ah as above. Let E be a given noninereasing, absolutely 
continuous function from [0, -|-oo) on [0, -|-oo) with E{0) > 0, M > 0 and j3 is a given 
parameter such that 


m 


</ 3 . 


2L(F'(r2)) 

In addition, E satisfies the following weighted nonlinear inequality 




2/3 


-)E{t)dt<ME{S), V0<5<r. 


(4.6) 


Then E satisfies the following estimate: 


> 


M 


0 


H'{rl 


(4.7) 


Furthermore, if limsupA//(x) < 1, then E satisfies the following simplified decay rate 

x—>0+ 


E{t) <2/3(77' 


1/kM\ 


t 




(4.8) 


for t sufficiently large, and where k > h is a constant independent of E{0). 
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5 Proof of main results 


In this section, we prove the main results including Theorems 13.Ill3^ and the decay rates 
in Examples 13.1113.41 


5.1 Proof of Theorem 13.11 

We define the following unbounded nonlinear operator ^ in by 

AU = {p, Au- aq - p{.,p),q, Av + ap ), 


with the domain 

D{A) = {U AU 

It is easy to check that D{A) = H (Hq) x Hq)^. Moreover, since p is nondecreasing with 
respect the second variable, we have for all U,U G D{A), 

{AU — AU ,U — U) = — f {p{x,p) — p{x,p)){p — ^dx < 0, 

Jn 

Thus —A is a monotone operator. We now claim that —A is a maximal operator. We 
proceed as follows. We denote by A the unbounded operator in defined hy A = —A and 
D{A) = UHq. Then / —^ is invertible as an operator acting from Hq in so that 

the operator {I—A)~^ is a well-defined, self-adjoint and if tc G then {I — A)~^w G 
Then for any F = (/, g,h,r) G "H, the equation 


{I-A)U = F 

with U = {u,p,v,q) G D{A) is equivalent to 

u — Au -|- a{I — A)~^{au) + p{.,u — f) = G, 
V = {I - A)~^{H 2 + au ), 

P = u- f ,q = v-h, 


where 


We define for 9 G 


Hi — 9 F f A ah G , 

H 2 = r F h — af G L? , 

G = Hi-a{I- A)-^H 2 G . 


f 

R{x,9) = / p{x,s)ds. 

Jo 

Let us define the functional J : Hq i—M defined by 

J{u) = J + |Vttp-|- {I — A)~^^^{au) ) + R{x,u — f) — Gu^dx . 


(5.1) 


(5.2) 
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Note that thanks to our hypotheses, \p{x, s)| < (^(l + |s|) for all (x, s) G 17 x M, so that J is 
well-defined and continuously differentiable on Hq. Moreover, we have 

J'{u).(p = f {uip + Vu ■ V(/? -|- (/ — (/ — + p{x, u — f)(p — Gip)dx . 

Jn 

We denote by 11 • 11 the norm. Since p is nondecreasing with respect to the second variable, 
J is a convex function and we also have 

so that J{u) —> -|-oo as ||Vtt|| —> -|-oo. Hence J is coercive. Therefore J attains a minimum 
at some point u G Hq, which satishes the Euler equation 

J'{u) = 0. 

The usual elliptic theory implies that the weak solution u of the variational problem 

iu^ Hi, 

= 0, ^ip£Hl 

is in . Hence u G H Hi. By dehning v as in (15. Hi , and p ,q as in (El]), it follows that 
U = {u,p,v,q) G D{A) and {I — A)U = F. Hence —A is a maximal monotone operator. We 
conclude Theorem 13. II using the classical theory of maximal monotone operator (see e.g. [21| 
and the references therein). 

5.2 Proof of Theorems 13.21 and 13.31 

The proof will be divided in three steps, following those described in m (see also mm- 

Step 1 : we hrst prove that the energy E satishes a suitable dominant energy estimate. 
This is the step in which the geometric assumptions PMGC on both the damping and 
the coupling regions are used, together with suitable multipliers adapted to the coupled 
structure of the wave-wave system. The proof is valid without specifying the growth 
assumptions on the feedback p. 

Step 2 : we then prove that nonnegative and nonincreasing functions E satisfying a suitable 
dominant energy estimate, satishes a general weighted nonlinear inequality. In the 
case of polynomially growing feedbacks p, the proof is easier since the weight function 
for integral inequalities is known. The general growing case relies on the optimality- 
convexity method of the hrst author [5|. 

Step 3 : we deduce energy decay rates, applying Corollary 14.11 for polynomially growing 
feedbacks, whereas applying Theorem 14.21 for general growing feedbacks. 
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Let us start with Step 1. We use the dominant energy method as developed and explained 
by the first author in UM- This method consists in estimating time integrals of the nonlinear 
weighted energy of the system by corresponding dominant weighted energies, here in the 
frictional case, it means by respectively the nonlinear kinetic energy and the localized linear 
kinetic energy. Note that this step is valid for feedbacks with polynomial as well as arbitrary 
growth close to the origin. 


Theorem 5.1. [Weighted dominant energy method] We assume that (IHCp holds where ujc 
satisfies the PMGC and that p € (7(0 x M) is nondecreasing with respect to the second 
variable, and p{-,0) = 0 on 0. Let uj^ be a given subset of Ll satisfying the PMGC. Let 
(j) : [0, +oo) (-)• [0, +oo) be a non-increasing and absolutely continuous function. Then, there 
exist constants 6i > 0{i = 1,2,3) and a* > 0 depending only on Ll,ujd but independent of 4> 
such that for any initial data in the energy space, for all a"*" G (0, a*], the total energy of the 
system (ini) satisfies the following nonlinear weighted estimate 

rT 

(f){t)E{t)dt < 5i(f){S)E[S) + 

rT r rT 


f 


^2 [ ([{t) f p^{x,u’)dxdt + f 4>{t) [ \u'\^dxdt, 

Js Jn Js Jujd 


(5.3) 


Before proving Theorem 15.11 we give a Lemma on a weighted energy estimate for a non- 
homogeneous wave equation. The proof of Lemma l5.II will be given in Section [6l 


Lemma 5.1. Letu be a nonempty open subset of LI satisfying the PMGC. Let : [0,+oo) e-)• 
[0, +oo) be a non-increasing and absolutely continuous function. Then, there exist con¬ 
stants rji > 0 {i = 1,2,3, 4) independent on (j) such that for all G Hq x L‘^ , all 

f G T^((0, oo); L^(0)) and all 0 < S <T, the solution u of 

u” — An = / in 0 X (0, +oo), 

< n = 0 on rx(0,+oo), (5-4:) 

(n,n')(0) = (n°,n^) in LI 


satisfies the estimate 


f 


(j){t)e{t)dt 


< r]i(j){S)[e{S)-\-e{T)]r ]2 -fi'{t)e{t)dt 

Js „ 


+73 


/ m [ I/p 

Js Jn 


dxdt J-r ]4 / 4>{t) / \u'\‘^dxdt. 


where e{f) := ^ Jf^(|n'p + \Vu\‘^)dx. 


(5.5) 


Proof of Theorem 15.11 

We first consider smooth initial data, then system dni) admits a unique solution (n, v) G 
C([0, Too), n Hlf) n iyi’-([0, +oo), {Hlf). 
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Let the weight function be a non-increasing absolutely continuous function. Let then 


1 


ei{t) ■■= - J^{\u I + |Vu| )dx. 

We now apply Lemma l5.II to the first equation of dni) with ui = ujd, f = —p{-, u') — a{x)v', 
with E given by (ll.2p and e{t) = ei{t). Using ei(t) < E{t) and the property that E is 
nonincreasing, we obtain for all 0 < S' < T and some constants rji > 0 (i = 1, 2,3,4) that 

f (f){t)ei{t)dt 

Js ^ 

< rii(l){S)[ei{S) + ei{T)] + Tj 2 [ {-^'{t))ei{t)dt 

Js _ 


+ 2 r ]3 f (l){t) f {x)\v'\'^dxdt + 2 r ]3 f f p'^{x,u')dxdt 

Js Jn Js Jn 


+rjA / ‘Jit) I \u'\‘^dxdt 


(5.6) 




< 


Cicj){S)EiS) + 2 ri^ f Jit) f p^ix,u)dxdt 

j, Js Jn ^ 

/ Jit) / \u'\‘^dxdt + 2 ps / Jit) / ix)\v'dxdt. 

Js Jtod Js Jn 


We now set 


1 


e2(i) := 2 I + I'^^l )dx , 


and secondly, we apply Lemma [541 to the second equation of dni) with uj = ujc, f = a(x)rt', 
with E given by ( 11 . 21 ) . e(t) = 62 (t). Again, using the inequality 62 (t) < Eit) and the property 
that E is nonincreasing, we obtain for all 0 < S < T and some constants 7 ^ > 0 (z = 1,2, 3,4), 

[ Jit)e2it)dt 

Js ^ 

< 'yiJiS)[e2iS) + e2iT)] + ^2 f i-j'it))e2it)dt 

fX' r ^ f ( 6 - 7 ) 

+73 / Jit) / cy^ix)\u\‘^dxdt + / Jit) / \v'\‘^dt 

Js Jn Js Jljc 

< C2(/>(S)U(S)-|-74 f Jit) f \v'\‘^dxdt +'Js f Jit) f a^ix)\u'\‘^dxdt. 

Js Juj^ Js Jn 

Let > 0 be a real parameter to be chosen later on. Adding (15.6D to 5- ( 1521 ), we obtain 
that for all 0 < S < T and all <5 > 0 

[ Jit)ieiit) + 6e2it))dt 

Js ^ 

< Cil + 6)JiS)EiS) +C [ Jit) [ p^ix,u')dxdt 

T . ' 6 f 

-|-C / Jit) / \u dxdt + Ciaj^-\ -) / Jit) / a\v'\^dxdt 

Js ^ Jwd O'- Js Jn 

+C5aj^ / Jit) / a\u'\^dxdt, 

Js Jn 
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where C denotes generic positive constants which may vary from one line to another. Next, 
we estimate the term fg 4>{t) J^a{x)\v'\‘^dxdt through the coupling relation. Obviously, the 
following identity holds for the solution (u, v) of system dni): 

[ 4>{t) [ [v'{u" — An + a{x)v' + p{x, n')) + n'(n" — An — a{x)u')]dxdt = 0. 

Js Jn 

After integration by parts, we obtain by Cauchy-Schwartz inequality that for all > 0 
[ (f){t) f a{x)\v'\‘^dxdt 


f (p{t) f a{x)\u'\‘^dxdt — 4>{t) f {u'v' + Vu-Vv)da 

Js ^ Jn i Jn 

+ / / (u'v'+ Vu ■ Vv)dxdt 

Js^ Jn 

~ / / p{x,u')v'dxdt 

Js Jn 


(5.9) 


< / (p{t) I a\u'\^dxdt + C(p{S)E{S) 

rT 


IS ^ Jn 

rT 


+ei 


^(t) f \v'\‘^dxdt + — [ (t>{t) f p‘^{x,u')dxdt. 

Jn Js Jn 


Using ()5.9p in (15.8h . we obtain for all 0 < 5 < 1 and all > 0 that 
[ (f){t)[ei{t) + 5e2it))dt 

Js ^ 

< C{l + 5 + aj^-\ - )(/)(5')£'(5) -|--|-(a+H-)— ^ J 4>{t) j p^{x,u')dxdt 

+C [ [ \u'\‘^dxdt + C 2 (a^ +—)ei f 4>{t) [ \v'\‘^dxdt 

Js Juid ^ Js Jn 

+Ciaj^{aj^-\ -) [ (pit) [ lu'l'^dxdt, 

a-' Js Jn 

where C, Ci and C 2 are generic positive constants. Thus, 

(1 - 2C'iQ;+(a+-I-—)) [ (p{t)ei{t)dt + [6 - 2C2ei(a++—)) [ (p{t)e 2 {t)dt 
«- Js «- Js 

r A 1 r'^ r 

< C'(1 + 5 + Q!+-|- )(/)(S')£'(5)-|-+ (a+-|-)—^ / <p{t) / p^{x,u')dxdt 

+C [ (p{t) [ \u'\ 

Js Ju)d 

Let q;+ be small so that 

We then fix (5 > 0 so that 

and choose next ei > 0 so that 


(5.10) 


‘dxdt. 


(5.11) 


0 < a+ < Cl := 


4Ci 


(1 - 2C'ia+(a+ —)) > - 


b - 2C2£i(a+ + —)) ^ - . 
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With these successive choices of a+,5 and ei, we deduce that 


r (j)it)E{t)dt < C{l + 5 + a+ + ^)cj){S)E{S) 

+C'fl+(a+H-)—) [ (pit) [ p^(x,u')dxdt (5.12) 

^ ^ o— Js Jn 

+C f (pit) f \u^\‘^dxdt. 

J S ^ d 


For initial data in the energy space, we conclude by a standard argument using density of 
n Hq) X Hq in Hq X L^, together with the dissipativity of the underlying nonlinear 
semigroup. This ends the proof of Theorem 15.11 □ 


Remark 5.1. The constant Ci, and therefore the constant a* depends only on H , 0 ;^. 


Proof of Theorem 13.21 [Case of polynomially growing feedbacks] 

Step 2. Assume that (HFp) holds. 

We first consider smooth initial data, then system admits a unique solution {u, v) G 
(^([O,+ 00 ), n Hq)^) n VF^’°°([0,+ 00 ), (Hq)^). Moreover, the total energy E defined by 
(na) is absolutely continuous and is non-increasing due to the monotonicity of p, that is 


E'{t) = — p{x,u')u'dx < 0. 


to. 


(5.13) 


Let t ^ 0 be fixed and .— ^x G . [?/ (t,x)| ^ 1} and ^ * [^ (^?^)l — 

1}. In short, we just write in the sequel. Then it follows from (HFp) and ()5.13p that 

[ (p{t) [ \u'\‘^dxdt <-C [ (p{t)E'{t)dt < C(p{S)E{S) (5.14) 

Js Juj° Js 


Similarly we obtain from (HFp), (I5.13P and Young inequality that for every 62 > 0, 


r (Pit) j i^'i 

J S J u)\ 


^'dxdt<C / (p{t) / \p{x^u)u\p-^^dxdt 


2 

p+i 


<(7/ ( I \p{x^u)u\dx 


<C (pit){-E'{t))^dt 

Js 

[ [e2{4iit))^ - C{e2)E'{t)]dt 

Js 

82 [ {(Pit))f^^dt + Cie2)EiS). 
Js 


dt 


(5.15) 


< 


< 


where 0 ( 82 ) > 0 stands for a constant depending on 82 (going to -|-oo as 82 goes to zero). 
Summing (I5.14p and (I5.15p gives 

f (pit) f \u'\‘^dxdt < max (ci 82 ),C 4 '{S)')EiS) + 82 [ i(pit))p-^dt, (5.16) 

Js Jlj^ ^ ^ Js 
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Similarly, let := {x £ : \u'\ ^ 1} and := {x € 17 : |ri'| < 1}. Note that these 

subsets depend as above on t. We get from (HFp), (|5.13l) and Young inequality that for every 
63 > 0 

r-T (■ rT 


[ (j){t) f p‘^{x,u')dxdt = [ ( [ p'^{x,u')dx+ [ p^{x,u')dx\ dt 

Js Jn Js \Jno Jq£ ) 

< C^{S')E{S') [ (f>{t) f \u p{x,u')\p+^dxdt 

Js Jn 

[cie3),C<PiS))E{S)+e3j^ m)f^dt. 


(5.17) 


< max 


We now choose the weight function (j) as follows 


-p-i 


= E—{t), Vt^O. 


(5.18) 


Combining (15.31) , (I5.16|) . (I5.17D . together with this choice for (j), and letting 62 , 63 small 
enough, we obtain that for all 0 < S' < T, 

rT 


f E'^{t)dt< max (Ci,C2E'^{0)'^E{S), 
J s 


where C'i,C '2 are positive constants independent of ill(O). 

We finish the proof by Step 3, which says that nonnegative, nonincreasing functions E 
satisfying a polynomial nonlinear integral inequality, then satisfy a polynomial decay rate. 
Applying Corollary 14.11 with a = > 0, and M = C{1 + E'^{Q)), we end the proof of 

Theorem 13.21 and obtain when i?(0) > 0 

E{t)<CEip)t-^/^^-^\ Vt ^ Ci.E^^(0) + C 2 , (5.19) 

where 


Ce{o) — 


(max(Ci,C2^'^(0))(l + ^))^^“. □ 


Proof of Theorem 13.31 [Case of general growing feedbacks] 

In Step 2 , We shall prove the following theorem 

Theorem 5.2. [Optimal-weight convexity method] Assume the hypotheses of Theorem 15. il 
Assume furthermore that [HFg) holds where g is such that the function H defined in (12. 3p is 
strictly convex on [OjTq], and g'{0) = 0. We define L by ( 12 . 4p . Let the initial data be in the 
energy space and be non vanishing, (u, v) be the solution of dni) and E be its energy. Then 
E satisfies the following nonlinear weighted integral inequality 

(tfA\E{t)dt<ME(S), V 0 <S<T, (5.20) 


f 


L 


-11 


where [3 and M are respectively given by 

/3 = max (^C 2 , 


m 


2L{H'{r^,)) 
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and 

M = 2Ci{l + H'{r^)), 

with Cl > 0,6*2 > 0 depending on 6i (i = 1,2,3), Q, ujd but independent of E{0). 


Remark 5.2. This method is called the optimal-weight convexity method according to the 
property that the weight function f is chosen in an optimal way by setting 




2/3 


thanks to suitable convexity arguments relying both on Jensen and Young’s inequalities for 
an appropriate convex function. 


Proof. We consider as before smooth initial data, then the solution (u, v) of (ini) is in 
C([0, +oo), n n +oo), (Rq)^). Moreover, the total energy E satisfies the 

dissipation relation (15.131) . Thanks to Theorem 15.11 we know that E satisfies the weighted 
dominant energy estimate (15.3p . We shall now use the optimal-weight convexity method of 
the first author [S] to build an optimal weight function f to prove that the two terms 

f (l){t) f p‘^{x,u')dxdt and f (^{t) f \u'\‘^dxdt 

Js Jn Js Jtod 

in (|5.3I) . are bounded above by the term 

CE{S){l + (f){S)) + C [ E{t)(j){t)dt V0<5<r. 

Js 

We proceed as in [5]. Choose a parameter eo sufficiently small, e.g. eo = g{rQ)). 

For fixed t ^ 0, we define the subset = {x G 12, \u'{t, x)| < eo}- Now thanks to (HF^), 
we know that (|2.2I) holds. Hence, since g is increasing, we have 

\p{x,u'{t,x))\. 


9{- 


K 


-) < \u\t,x) \ , for a.e x G 12} , 


(5.21) 


where K = C2||a|loo) with || • l|oo standing for the norm. Now, we can note that parameter 
eo bas been chosen to guarantee the following two properties 

1 f |p(x,u'(t,x))p 


-dx G [0,ro], 


and 


1^11 ln\ 

Wk L, p{x, u {t, x))u'{t, x)dx G [0, R(xo)], 


(5.22) 


(5.23) 


hold. Since H has been assumed to be convex on [0,rQ] and thanks to (I5.22p . the Jensen’s 
inequality, and (I5.21|) . we obtain 


H 


^ 1^11 Iq\ 


\p{x,u'(t,x))\‘^ 


iF2 


""'-Mile 


\p{x,u'{t,x))\‘^ 

iF2 


dx = 


\Q\\K 




\p{u{t, x))| g{ )dx < {t,x)p{x,u {t,x))dx . 
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By (I5.23j) . we deduce that 

rT 


[ (j){t) f \p{x,u'{t,x))\'^dxdt < f K'^\Q,\\(p{t)H ^ f , [ u'{t,x)p{x,u\t,x))dx')dt, 

Js Jn\ Js Jn{ ^ 

and using further Young’s inequality, the dissipation relation (j5.13j> . we obtain 


r-T /> i-T 


/ \p{x,u'{t,x))ydxdt<K^\n\ H*{4>{t))dt + KE{S), V0<5<T. 

's Jq\ Js 


<T. 


On the other hand, we prove easily as in the proof of Theorem 13.21 that 

f 4>{t) f \p{x,u {t,x))\‘^dxdt < KE{S)(j){S), V 0 < S' 

Js Jn\n\ 

Adding these two inequalities, we obtain 

[ [ \p{x,u'it,x))\‘^dxdt<K^\n\ [ H*{(J){t))dt + KE{S){l + ^{S)), V0<S<r. 

Js Jn Js 

(5.24) 

We now turn to the estimate of the localized weighted linear kinetic energy. Thanks to (HFg ), 
we know that ()2.2I) holds. Choose a parameter ei sufficiently small, e.g. ei = min{ro, ^(ri)} 
where ri is defined by 

rl = (;|^(^o)) > 

where k = cia_. For fixed t ^ 0, we define the subset = {x € , \u'{t, x)\ < ei}. Thanks 


to ()2.2p and since (HFg) holds, we have 

g{\u'{t,x)\) < 


\p{x,u' {t,x))\ ^ 

, for a.e x £ uE . 
k 


(5.25) 


Now, we can note that parameter ei has been chosen to guarantee the following two properties 

1 




\u'{t, x)\^dx e [0, Tq] , 


and 




p{x,u'{t,x))u'{t,x)dx G [0, Ff(rg)] , 


(5.26) 


(5.27) 


hold. Since H has been assumed to be convex on [0,rQ] and thanks to (|5.26p . the Jensen’s 
inequality, and (I5.25p . we obtain 


H 


■ — f- [ |u'(t,x)p(ix') <- — f- [ H{\u'{t,x)\'^)dx 

K\ Ji^i ’ K\ Jooi 


UJ 


\u'{t, x)\g{\u'{t, x)\)dx 


d\ 


< 


- — f— [ u'{t,x)p{x,u'{t,x))da 


Thanks to (15.271) . we deduce that 


[ [ \u (t,x)\‘^dxdt < [ ^(- —7— [ u'{t,x)p(x,u'(t,x))dx')dt. 

Js Js ^\^d\>^Jn\ > 
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and using further Young’s inequality, the dissipation relation we obtain 

[ [ \u\t,x)\‘^dxdt < \ujd\ [ H*{(l){t))dt +—E{S), V0<5 

Js Js k 

On the other hand, we prove easily as in the proof of Theorem 13.21 that 


<T. 


r m [ 

Js Ju)d\uj\^ 


1 


u'{t,x)\^dxdt < —E{S)4>{S), V 0 < S' < T. 


k 


Adding these two inequalities, we obtain 


[ (l){t) [ \u'{t,x)\'^dxdt < \ujd\ [ H*{(j){t))dt->r—E{S){1 + (I){S)), V 0 < S < T . (5.28) 

Js Juid Js k 

Using (I5.24P and (15.281) in the weighted dominant energy estimate (15.31) . we obtain 


L 


(J){t)E{t)dt < CiE{S)il + (j){S)) + C 2 [ H*{(p{t))dt, V 0 < S < T, (5.29) 

IS Js 

where the constants Ci,C 2 depend only on the 6 i for i = 1,2,3 and on |fl| and \iOd\ in an 

explicit way. In particular, they do not depend on cj). 

Let 


13 = max (^C 2 , 


m 


where L is defined in (j2.4l) . Since Li is a nonincreasing function, and thanks to (12.5p . we have 

E{t) ^ E{0) 


< 




2/3 - 2p 

Hence, since L~^ is defined from [0,rQ) onto [0,+oo), we can define 4> by 


<ti{t) = L > V t ^ 0. 


(5.31) 


By definition of L, 0 is a nonnegative, non increasing and absolutely continuous function on 
[0,+oo). We hrst note that 

(/{S) <H'{rl), VS^O. (5.32) 

Then, thanks to our ’’optimal” choice of the weight function </) and to the definition of L, we 
have 




2/3 




This implies 


c2H*m)) <pH*m)) = vt^o. 


Combining this estimate together with (|5.32l) in (I5.29p . we obtain 

(J^'\E{t)dt<ME{S), V0<S 


f 


L 


-1 / 


<r. 


where 


M = 2Ci{l + H'{rQ). 


(5.33) 

(5.34) 

□ 
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We finish the proof of Theorem 1,3.31 bv Step 3. Thanks to Step 2, E satisfies the weighted 
nonlinear integral inequality (I5.33p . where /3 is dehned by (|5.30p . M is defined by (I5.34p . 
Hence applying Theorem 14.21 with this /3 and M we deduce that E satishes the decay rate 
(14.71) in the general case, and the simplihed decay rate (14.81) if limsupA//(x) < 1. This 
concludes the proof of Theorem 13.31 □ 

5.3 Proof of Theorem 13.41 

Thanks to our hypotheses, the simplified upper energy estimate (|4.8p of Theorem 14.21 
holds, so that E{t) converges to 0 as t goes to infinity. Hence, there exists Tq ^ 0 such that 

<(-)', Vt>ro, (5.35) 

where % = 4y^i?i(0). Hence 

-is\/E{t) G [0, rg] for all t ^ Tq . (5.36) 

On the other hand, thanks to the regularity of u (see [TO] for details), we have 

lk'(i, OIIl'-co) < ls\/E{t) 

Using this inequality in the dissipation relation 

—E'{t)= / a{x)u'g{u')dx , 

Jn 

together with (I5.36D . we deduce that for all t ^ Tq 

-E'{t) < — E{t)H{-fs\/E{t)). 

Is 

Therefore, we have 

E{t)>(^^K-\\\a\\Loo^n){t-To)))\ yt^To, (5.37) 

where K~^ denotes the inverse function of K defined by 

= VrG(0,v^:^), (5.38) 

where 

zo = 'ysVE{To). (5.39) 

We denote by z the solution of the following ordinary differential equation 

z{t) + \\a\\L^{n)H{z{t)) = 0,2;(0) = zo . (5.40) 

Then we have the relation 

z{t-To) = K-\\\a\\Lo.^n){t-To)), ^ t ^ Tq . (5.41) 

We now use the following comparison Lemma, that we recall for the sake of completeness. 
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Lemma 5.2. Lemma 2-4] Let H be a given strictly convex function from [OjTq] to M 
such that H{0) = H'{0) = 0, where r^ > 0 is sufficiently small and define Ah on (OjTq] by 

Let z be the solution of the ordinary differential equation: 

z'{t) + KH{z{t)) = 0 , z{0) = zq t>0, (5.42) 

where zq > 0 and k > 0 are given. Then z{t) is defined for every t > 0 and decays to 0 at 
infinity. Moreover assume that (HFi) holds. Then there exists Ti > 0 such that for all R > 0 
there exists a constant C > 0 such that 

z{t)>C{HT^(^j), yt>Ti, (5.43) 

where Ti is a positive constant. 

We apply this Lemma to the solution z of (15.401) with R = 1 and k = ||a||ioo(Q). Thus, 
there exist two constants Ti > 0 and > 0 such that 

z{t)>Cs{H')-^(^^'^ , yt>Ti, (5.44) 

Combining (j5.44p together with (I5.37p and (j5.4ip . we obtain the lower estimate (13.81) . □ 

Remark 5.3. The constant C of the above Lemma 15.21 depends explicitly on k, R (and in 
addition of p, if the second alternative of {HFi) holds). This dependence is given in the proof 
of Lemma 2.4 in [8]. Moreover, one may assume that tq = oo in Lemma 15.21 In this case the 
interval [0, Tq] becomes [0,+oo). 

5.4 Proof of the decay rates given in Examples 13.1113.41 

Example 13.11 We have H{x) = for x G [0, Tq]. Thus H'{x) = and H 

is strictly convex on a right neighborhood of 0. Moreover, Ah{x) = < 1 for all x € [0, Tq]. 

We easily conclude applying ()4.8p of Theorem 14.21 and Theorem 13.41 for the lower estimate in 
the one-dimensional case. 

Example 13.21 We have H{x) = ^Jxe~ for x G [0,Tq]. Thus H'{x) = ^^ {\ + i), and 

H is strictly convex on a right neighborhood of 0. Moreover, we have Ah{x) = i ^ i for all 

___ ' 2 ) 

X > 0 sufficiently close to 0, so that lim Ah(x) = 0. We apply (14.8p of Theorem 14.21 So we 

X-I-0+ 

set x{t) = (R')”^^^^y Then one can prove that x{t) is equivalent to in(f)-in(KM) ^ Soes 
to -boo. We therefore obtain the desired upper bound, using this equivalence. One can show 
that the second alternative of (HFi) holds for any p > 1 (see subsection 7.10 in [S]). Thus, 
we obtain in the same way by Theorem 13.41 the lower estimates in the one-dimensional case. 

Example 13.31 We have F[{x) = ^/x)Y for x G [0,rQ]. Thus H'{x) = 

i^(p-i)/2(in(i/,^^))ij g(ln(l/y^))“^^ and H is strictly convex on a right neighborhood 

of 0. Moreover, Ah{x) = for all x > 0 sufficiently close to 0, so that 
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lim Ah(x) = We apply (14.81) of Theorem 14.21 So we set x(t) = iH') and 

^>0+ \ ^ / 


yi't) = • Then one can prove that is equivalent to 

as t goes to +oo. On the other hand, computing ln(x(t)) and ln(y(t)), we find that ln(x(t)) 
is equivalent to ln(y(t)) as t goes to +oo. Using this relation in the previous one, we find 
that x{t) is equivalent to Dt~'^/^P~^^ ("in'"''' 


(p-l)/2 , 


P+1 


, where D is an explicit positive constant 
which depends on k, M,p and q. We therefore obtain the desired upper estimate, using this 
equivalence. We obtain by Theorem 13.41 the lower estimates in the one-dimensional case. 
Example 13.41 We have H{x) = for x £ [0,rg]- Thus, we have H'{x) = 

-|-p^ln and H is strictly convex on a right neighborhood 

of 0. Moreover, we have Ah{x) = - - - - -for all x > 0 sufficiently close to 0, so 

l+p(ln(l/Vx)j 

that lim Ah(x) = 0. We apply (14.8p of Theorem 14.21 So we set x(t) = ^'] and 

3:^0+ \ r / 


y{t) = Then one can prove that ln(x(t)) is equivalent to ln(y(t)) as t goes to 

-|-oo. We further set z{t) = ln(l/yx(t)) so that z{t) goes to +oo as t goes to -|-oo, then 
we have zP{t){l — 9{t)) = ln(2^), where 6{t) = z^~P{t) + ln(l + pzP~^{t))z~P{t), so that 
6{t) goes to 0 as t goes to +oo. Hence we have x{t) = e (i-0(t))i/p _ ’^g check 

that ln(2^)^^^(l — (1 — 9{t)y/P) goes to 0 as t goes to +oo. Hence, x{t) is equivalent to 
g-2(in(t))i/p g^g ^ goes to +00. We therefore obtain the desired upper estimate. One can show 
that the second alternative of {HFi) holds for any p > 1 (see subsection 7.10 in [H]). Thus, 
we obtain by Theorem 13.41 the lower estimates in the one-dimensional case. 


6 Proof of Lemma 15.1 


In this section, we prove Lemma 15.1 1 bv the piecewise multiplier method which relies on the 
geometric assumptions PMGC on the subset cu C H. Denoting by and Xj (j = 1, • • • , J) 
the sets and the points given by PMGC, we have uj D Ni.{i_}'j^pyj{xj) U (H \ n H. 

Here, N^{U) = {x £ M"', d{x, U) < e} with d(-, U) the usual euclidian distance to the subset 
U of M”, and 7j(xj) = {x £ Pj, (x — xj) ■ Uj > 0}, where Uj denotes the outward unit normal 
of the boundary P^ = dQ,j. 

Let 0 < £o < £i < €2 < £ and define Qi := N^^[{J’j^^'yj{xj) U (H \ = 0,1,2). 

Since (Hj\Qi) n Qo = 0) we introduce a cut-off function ipj £ satisfying 

0 < V'j ^ 1) V'j = 1 on Qj\Qi; ^|)j = 0 on Qq. (6.1) 


For mj{x) = x — Xj, we define the vector field on H: 

'0j(x)mj(x) if X £ 0,j,j = 1, • • • , J 
0 if X £ H\ 


h{x) = 


Using the multiplier (j)(t)h{x) • Vu to equation (15.4p : 


/ Hi) [ hix) ■ V«( 

Js JQj 


u" — Au — f)dxdt = 0 


( 6 . 2 ) 
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leads to 


/ 1 

—h ■Vu+-{h- z/,)(|u'|2 - \Vu\'^)dTdt 
r, duj 2 

-\T rr. 


’ 4>'{t) / u'h-Vudxdt 

S Jui 


Is 

4){t) / u'h ■ Vudx 

J 

- f (pit) f -divhilu'l"^ — \Vu\‘^)dxdt 

Js Jui 2 

- r<pit) [ 

Js Jn 


(6.3) 


E du du dhk , ^ r 

dxidxkdxi ■ 


dxdt. 


Thanks to the choice of ipj, the terms in the left hand side of (I6.3h vanish except on the 
boundary (Tj \ 7 j(xj)) n T. Since n = 0 on this part of boundary, then u' = 0, Vn = ^ Vj. 
Hence, the left side of (16.3p becomes 

<'T r 2 


Us 




dus 


ipjfrij ■ KjdVdt < 0. 


Therefore (16.311 implies that 


< - 


f (pit) f -divkilu^l"^ — \Vu\‘^)dxdt + f (pit) f 
Js Juj 2 Js Jui 

-1 T ^ 

(pit) f u'h-Vudx + f (p'it) f u' 

J Q j J S J 

r m [ 

Js Jn, 


‘3 i,k 
h ■ Vudxdt 


du du dhk 
dxi dxk dxi 


dxdt 


h ■ Vtt fdxdt. 


Since, moreover /i(x) = rrijix) on we obtain that 


(pit) 


N 


ifljXQi 


— (|ri'l — |Vuf) + \Vufdxdt 


< - 


+ 


4>it) [ u' 
Jfij 

f m [ 

Js Jq 

fu 


h ■ Vudx 


+ 


r^it) [ n' 

Js Jfli 


h ■ Vudxdt 


h -Vu fdxdt 
1 


Qj nQi 


divhi\u\‘^ — |Vnp) + ^ 


2 ,fc 


du du dhk 
dxi dxk dxi 


dxdt. 


Summing the above inequality on j and using the facts that H\(5i = and 

hix) = 0 on H \ we obtain 


< - 


[ 4’i't) [ ^iW\‘^ — |Vn|^) + \Vu\‘^dxdt 

Js Jn\Qi 2 

(pit) / u'h ■ Vudx 

rp 

I (pit) / h-Vu fdxdt 
s Jn 


+ f (p'it) f u'h-Vudxdt 
^s Js Jn 


(6.4) 


rm [ 

Js JnnQi 


hf\ '|2 1 ^ |2n ^ du du dhk 

-d,vk(]u\ -|V»| ) + ^___ 


i,k 


dxdt , 
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Using the second multiplier —^4>{t)u for (15.4^ : 


N -1 


[ [ «( 

Js Jn 


yields that 


A^- 1 


+ 


2 

N - 1 
2 

(N-l) 


u" — Au — f)dxdt = 0 . 


f f “ \u'\‘^)dxdt 

Js Jn 


4>{t) / u'udx 


rT 


(6.5) 


+- 


2 

N - 1 




4>'(t) / u'udxdt 


/ (/)(t) / u fdxdt. 

Js Jn 


Adding (16.5D to (16.41) and using Cauchy-Schwarz and Poincare’s inequalities, we obtain for 
all (5i > 0 that 


(j){t)e{t)dt 


|Vrip + \u‘ 


dxdt 
T 


= r m 

Js ■ ^ 

< C(l){S)[e{S) + e{T)]+C I -<p'{t)e{t)dt 

T H 

+(5i f (j){t)e{t)dt + Y f </’(0 f Ifl"^dxdt 

J ^ J S J ^ 

+c [ <^(t) 

Js 


( 6 . 6 ) 


(|u'p + \Vu\‘^)dxdt 

's JnnQi 

Compared to the desired estimate (15.51) . the term concerning |VmP on the right hand of 
(j6.6l) is crucial. We just follow the techniques developed in [27| to deal with this term. 

Since E.^\Q 2 CiQi = ij), there exists a cut-off function ^ e (^(^“(M) such that 

0<^<1, C = 1 on Qi, ^ = 0 on M^\Q 2 . (6.7) 

Applying now the multiplier (j)(t)^{x)u to (15.4p gives, after integration by parts, that 

[ 4^ij) f ^\'^u\‘^dxdt 

Js„ Jn 

1 


f J’i't) I + ^u^A()dxdt 

Js Jn J 

j-T 

+ 


f 4>'{t) f ^uu'dxdt — (j){t) f ^uu'dx 

Js Jn V Jn . 


+ 


f 4>{t) f fdxdt. 

Js Jn 
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Then it follows from the dehnition of ^ that 


f (p{t) f \Vu\‘^dxdt < f (p{t) f ^\Vu\‘^dxdt 

Js JnnQi Js Jn 


< Cct){S)[e{S) + e{T)]+C -(l)'{t)e{t)dt 

Js ( 6 . 8 ) 

^dxdt 


+c [ m [ I/p 

Js^ Jn 

+C f (p{t) f + \u\‘^)dxdt. 

Js JnnQ 2 


Now it remains to estimate the term concerning jup in (j6.8p . Since R^\ujr\Q 2 = 0, there 
exists a function /3 G C'q°(M) such that 


0</3<l, /3 = 1 on Q 2 , /3 = 0 on . (6.9) 

Fix the t variable and consider the solution z of the following elliptic problem in space: 

{ Az = (3{x)u, in 17, 
z = 0 , on r. 


Hence, 2 and z' satisfy the following estimates 


|-2||l,2 < C 


f P{x)\u\‘^dx. 

Jn 

f fJ{x)\u'\'^dx. 

Jn 


\W\\h<c 

Applying the multiplier (j){t)z to (15.4p gives, after integration by parts, that 

rT 


( 6 . 10 ) 

( 6 . 11 ) 


f 4>{t) f /3{x)\u\‘^dxdt 

Js Jn ^ ^ 

4 >{t) / zu'dx — / / zu' 

^Jn 1 s Js Jn 

+ f (j){t) f [—z'u—zf)dxdt. 

Js Jn 


Hence, using the estimates (|6.1UI) - (l6.1ip in the above relation, and noting the dehnition of /3, 
we obtain for all (52 > 0 

cT r rT 


f <p{t) f \u\'^dxdt < [ <p{t) f f3{x)\u\‘^dxdt 

Js JnnQ2 Js Jn ^ 

< Ccl){S)[e{S)+e{T)] + C [ -(j)'{t)e{t)dt 

T ' ^ 

+<52 / (j){t)e{t)dt + — / (pit) / \u'\‘^dxdt 

Js ^ Js Jui 

t>2 Js Jn 


( 6 . 12 ) 


+- 


Inserting (16.81) and (I6.12p in (16.6p , then choosing hnally and 62 sufficiently small, we obtain 
(j5.5l) which ends the proof of Lemma l5.11 □ 
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